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Heterogeneous  Sorption  Processes  in  Subsurface 
Systems:  Deveiopment  and  Appiication  of  Modeis  to 
Batch  Systems 


Josqjh  A.  Pedit  and  Cass  T.  Miller 

Introdaction 

Contamination  of  sediment,  soil,  and  subsur&ce 
environments  by  hydro{diobic  organic  compounds  has 
led  to  significant  research  in  sorption  processes.  Many 
matfaematicat  models  have  been  formulated  for  model¬ 
ing  the  sorption  process;  most  use  a  deterministic 
fiamewmic  and  assume  an  idealized  homogeneous 
sorbent  structure.  However,  natural  systems  display 
heterogeneities  at  scales  ranging  fiom  the  megascale 
used  to  describe  large  natural  systems  to  the  microscale 
of  an  individual  grain  or  aggregate. 

Eiqperimental  investigations  of  sorption  processes 
are  often  conducted  in  batch  systems  in  an  efibrt  to 
isolate  the  sorption  process  fiom  advective  and  disper¬ 
sive  phenomena.  Modeling  of  such  systems  requires  a 
description  of  sorption  equilibrium  and  rate  properties. 
Most  sorption  models  assume  the  solid  phase  of  a  bulk 
sample  to  be  homogeneous  with  respect  to  sorption 
equililnium  and  rate  properties,  although  many  experi¬ 
mental  investigations  have  demonstrated  otherwise. 

This  article  provides  an  overview  of  models 
formulated  to  address  sorption  equilibrium  and  rate 
{xoperty  heterogeneities  within  a  bulk  sample.  A  subset 
of  tile  models  is  fit  to  an  experimental  data  set. 

Sorption  Modeb 

The  sorption  [xocess  is  assumed  to  be  difiiision 
limited.  The  time  required  for  a  solute-sorbent  system 
to  reach  equilibrium  has  been  observed  to  be  related  to 
the  magnitude  of  tiie  equilibrium  distribution  relation¬ 
ship  (Karickhoff  and  Morris,  1985;  Brusseau  and  Rao, 
1989),  and,  to  rome  extent,  particle  size  (Ball  and 
Roberts,  1991a,  1991b);  both  observations  are  consistent 
witii  a  difiusion-iimited  process.  Given  the  variations  in 
equilibrium  distribution  relationships  and  particle  size 
tiuti  are  possible  within  a  bulk  sample,  it  is  not  surpris¬ 
ing  that  sorption  modeling  apinxiaches  that  incorporate 
equilitnium  distribution  relationships  that  are  measured 
on  a  bulk  sanqile  and  incorporate  a  first-order  mass 


transfer  rate  model  have  met  with  only  limited  success 
in  batch  systems  (e.g..  Miller  and  Weber,  1986;  Wu  and 
Gschwend,  1986).  An  instantaneous  equilibrium 
fiaction  is  often  incorporated  into  rate  modeling  of 
batch  systems  to  improve  model  fits  (e.g..  Miller  and 
Weber,  1986;  Ball  and  Roberts,  1991b).  The  need  to  use 
an  instantaneous  equilibrium  fiaction  in  rate  modeling 
of  batch  systems  is  likely  to  be  a  result  of  model 
misspecification  caused  by  variations  in  sorption 
equilibrium  and  rate  properties  at  the  grain  scale. 

The  first  type  of  models  that  we  have  formulated  is 
an  extension  of  difiusion-type  models.  The  diffusion 
process  can  be  modeled  by  Pick’s  law  for  a  well-defined 
geometry  or  by  a  first-order  mass  transfer  ^proxima- 
tion. 

Diffusion  models  based  on  Pick’s  law  often  assume 
that  the  sorption  process  is  controlled  by  pore  or  surface 
diffusion  into  a  well-defined  particle  geometry,  usually  a 
sphere.  Pore  diffusion  is  envisioned  to  occur  in  the 
intrq)aiticle  pore  spaces  of  minerals  and  aggregates. 
Surface  diffusion  is  envisioned  to  occur  along 
intraparticle  pore  walls.  The  pore  diffusion  modeling 
^^nroach  has  been  extended  to  include  an  instantaneous 
equilibrium  fi»:tion  in  an  effort  to  improve  model  fits 
(e.g..  Ball  and  Roberts,  1991b). 

Our  general  model  based  on  Pick’s  law  incorporates 
the  major  features  that  have  been  used  by  other  re¬ 
searchers.  The  general  model  formulation  allows  for 
both  pore  and  surface  diffusion,  multiple  particle  classes 
with  different  physical  properties  and  sorption  param¬ 
eters,  and  an  instantaneous  equilibrium  fraction  for  each 
particle  class. 

Pirst-order  mass  transfer  model  formulations  for  the 
sorption  process  are  frequently  used  for  modeling  batch 
systems  (e.g..  Miller  and  Weber,  1986;  Wu  and 
Gschwend,  1986).  Wu  and  Gschwend  (1988)  attempted 
to  account  for  the  possible  effect  of  particle  size  varia¬ 
tions  within  a  bulk  sample  by  allowing  for  multiple 
particle  classes,  with  the  first-order  mass  transfer 
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Figure  1 .  Experimental  data  and  Toth  model  fits  for  the  distribu¬ 
tion  of  diuron  at  1 2  weeks  on  pulverized  and  unaltered 
Wagner  material. 

coefficient  for  each  particle  class  being  inversely  related 
to  the  square  of  the  particle  radius. 

A  fraction  of  equilibrium-Qqie  sorption  sites  are 
often  included  in  first-order  mass  transfer  models  wfrich 
have  been  used  to  model  batch  systems  (e.g.,  Karickhoff 
and  Morris,  1985;  Miller  and  Weber,  1986).  First-order 
mass  transfer  models  that  include  a  fraction  of  equilib¬ 
rium-type  sorption  sites  are  often  referred  to  as 
bicontinuum  models.  Physical  process  interpretations 
of  bicontinuum  modek,  often  referred  to  as  two-region 
or  mobile-immobile  fluid-phase  models,  assume  that  all 
of  the  sorption  sites  are  at  equilibrium  with  the  adjacent 
fluid-i^iase  solute  concentration.  However,  transfer  of 
solute  between  die  mobile  and  immobile  fluid  phases  is 
diffusion  controlled.  In  our  model  formulation,  the 
mobile  and  immobile  fluid  phases  are  viewed  as  the 
bulk  and  the  intr^iarticle  fluid  phases,  respectively.  A 
general  first-order  mass  transfer  model,  ^ch  allows 
for  multiple  particle  classes  with  equilibrium  and 
diftiision-cmitroUed  sorption  sites,  was  formulated  to  be 
consistent  with  the  general  model  based  on  Pick’s  law. 

The  two  general  model  formulations  outlined  above 
attempt  to  addre^  the  variability  that  is  possible  in  a 
bulk  sample.  E«;h  particle  class  of  a  bulk  sample  could 
be  assigned  (riiysical  properties  and  sorption  parameters 


with  these  model  formulations,  but  such  an  rqiproach  is 
impractical  from  an  experimental  and  modeling 
perspective.  A  more  viable  approach  is  to  describe  the 
sorption  equilibrium  and  rate  parameters  of  the  bulk 
sample  widi  a  continuous  bivariate  probability  density 
function,  which  we  have  tqiplied  to  a  first-order  mass 
transfer  framework.  In  our  model  formulation,  we 
assume  that  the  sorption  equilibrium  and  first-order 
mass  transfer  coefficients  are  random  variables  with 
their  distribution  given  by  the  bivariate  lognormal 
probability  density  function.  The  model  can  be  simpli¬ 
fied  when  the  sorption  equilibrium  parameter  is 
constant  throughout  the  bulk  sample.  The  simplified 
model  describes  the  continuously  distributed  first-order 
mass  transfer  coefficient  by  the  lognormal  probability 
density  function. 
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Figui«2.  Experimental  data  and  two-parameter  model  fits  for 
diuron  sorption  on  unaltered  Wagner  material. 


Application  of  Models 

A  subset  of  the  models  described  above  was  fit  to 
the  results  of  a  sorption  rate  experiment.  Sorption 
equilibrium  and  rate  eiqieriments  were  conducted  with 
the  herbicide  diuron  on  a  glacially-deposited  material 
(Wagner)  with  an  organic  carbon  content  of  0.22  ± 
0.02%.  Glacially-deposited  materials  are  typically 
heterogeneous  vrith  respect  to  particle  size  and  origin, 
and  the  Wagner  material  sample  is  no  exception. 

Equilibrium  eiqieriments  were  conducted  on 
unaltered  and  pulverized  material.  The  material  was 
pulverized  because  the  use  of  pulverized  material  in 
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soipti<Hi  equililMium  eiqperiinents  has  been  shown  to 
substantially  reduce  the  time  required  to  reach  equilib¬ 
rium  without  significantly  altering  the  final  equilibrium 
(Ball  and  Roberts,  1991a,  1991b).  The  experiments 
were  allowed  to  equilibrate  for  12  weeks.  The  results  of 
die  equilibrium  experiments  on  the  unaltered  and 
pulverized  materials  along  with  the  Toth  equilibrium 
model  fits  are  shown  in  Figure  1 .  The  pulverized 
material  sorbed  about  2.3  times  as  much  as  the  unal¬ 
tered  material  at  a  given  equilibrium  bulk  fluid-phase 
solute  concentration.  These  results  suggest  that  the 
unaltered  material  was  far  from  equilibrium  even  after 
12  weeks.  The  equilibrium  data  of  the  pulverized 
material  are  probably  more  representative  of  true 
equilibrium  and  the  Toth  model  fit  of  the  results  on  the 
pulverized  material  was  used  for  rate  modeling. 

The  rate  study  used  for  the  sample  application  of 
the  rate  models  h^  an  initial  bulk  fluid-phase  solute 
concentration  of  3050  mg^L*',  a  mass  of  solids  of  5  g, 
and  a  volume  of  bulk  fluid  of  24.9  mL.  Three  rate 
models  with  two  adjustable  parameters  were  fit  to  the 
data  set.  These  mo^ls  included  the  Fickian  diffusion 
model  with  a  single  particle  class  assuming  pore 
diffusion  with  an  instantaneous  equilibrium  fraction,  the 
first-order  mass  transfer  model  with  a  single  particle 
class  with  an  instantaneous  equilibrium  fraction,  and 
the  stochastic  version  of  the  first-order  mass  transfer 
model  with  a  continuously  distributed  first-order  mass 
transfer  coefficient  given  by  a  lognormal  fnobability 
density  fimction.  The  results  of  the  sorption  rate 
erqreriment  and  the  three  model  fits  are  shown  in  Figure 
2.  The  stochastic  version  of  the  first-order  mass  transfer 
model  provided  the  best  fit  to  the  e}q)eTimental  data. 

The  pore  diffusion  model  fit  was  only  slightly  worse,  but 
still  accurately  cultured  the  dynamics  of  the  system. 

The  single  particle  class,  first-order  mass  transfer  model 
with  an  instantaneous  equilibrium  fiaction  provided  a 
poor  fit  to  the  data.  This  model  overpredicts  early 
iqjtake  and  also  iq>proaches  equilibrium  too  rsq)idly. 

Obviously,  insight  into  the  proper  modeling 
qiproach  requires  additional  ejqterimental  information 
on  variations  in  sorption  equilibrium  and  rate  properties 
within  a  bulk  sample.  The  time  and  effort  involved  in 
such  detailed  experimental  work  would  be  prohibitive 
for  most  investigations  of  sorption  behavior.  Barring 
such  detailed  investigations,  modeling  sorption  rates  of 
bulk  solids  in  batch  systems  using  ai^noaches  similar  to 
the  single  particle  class,  pore  diffusion  model  with  an 


instantaneous  equilibrium  fiaction  or  approaches 
similar  to  the  lognormally  distributed  first-order  mass 
transfer  model  are  expected  to  (novide  good  descriptions 
of  experimental  results.  ^ 
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